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Abstract 

By modeling a linear polarizable and magnetizable medium (magneto- 
dielectric) with two quantum fields, namely E and M, electromagnetic 
field is quantized in such a medium consistently and systematically. A 
Hamiltonian is proposed from which, using the Heisenberg equations, 
Maxwell and constitutive equations of the medium are obtained. For a 
homogeneous medium, the equation of motion of the quantum vector 
potential, A, is derived and solved analytically. Two coupling func- 
tions which describe the electromagnetic properties of the medium are 
introduced. Four examples are considered showing the features and 
the applicability of the model to both absorptive and nonabsorptive 
magneto-dielectrics . 
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1 Introduction 

In a homogeneous and nondispersive medium, the photon is associated with 
only the transverse part of the electromagnetic field. In contrast in an in- 
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homogeneous nondispersive medium, the transverse and the longitudinal de- 
grees of freedom are coupled. In this case the quantization of the electro- 
magnetic field can be accomplished by employing a generalized gauge that is, 
V • (e(f)A) = 0, where e(f), is the space dependent dielectric function [1,2]. 
Using the gauge Y%j=i ~£r( £ ij(f)Aj) = 0, the generalization of this quanti- 
zation to the case of an anisotropic nondispersive medium is straightforward 
[3]. 

The quantization in a dispersive and absorptive dielectric represents one 
of the most interesting problems in quantum optics, since it gives a rigorous 
test of our understanding of the interaction of light with matter. The dis- 
sipative nature of a medium is an immediate consequence of its dispersive 
character and vice versa according to the Kramers- Kronig relations. This 
means that the validity of the electromagnetic field quantization in a nondis- 
sipative but dispersive media, is restricted to a range of frequencies for which 
the imaginary part of the dielectric function is negligible. Otherwise, there 
will be inconsistencies in electromagnetic field quantization process. 

In the scheme of Lenac for dispersive and nonabsorptive dielectric media, 
by starting from fundamental equations of motion, the medium is described 
by a dielectric function e(r, ui) without any restriction on its spatial behavior 
[4]. In this scheme, it is assumed that there are no losses in the system, so the 
dielectric function is real for the whole space. The procedure is based on an 
expansion of the total field in terms of the coupled eigen- modes, orthogonality 
relations are derived and equal-time commutation relations are discussed. 

Huttner and Barnett have presented a canonical quantization for electro- 
magnetic field inside an absorptive dielectric [5]. In their model, the medium 
is represented by a collection of interacting matter fields and the absorptive 
character of the medium is described by interaction of the matter fields with 
a reservoir containing a continuum of Klein-Gordon fields. In this model, 
eigen-operators for the coupled systems are calculated and electromagnetic 
field is expressed in terms of these operators. Also the dielectric function is 
derived and is shown to satisfy the Kramers-Kronig relations. 

Gruner and Welsch presented a quantization method of the radiation 
field inside a dispersive and absorptive linear dielectric starting from the 
phenomenological Maxwell equations, where the properties of the dielectric 
are described by a permitivity consistent with the Kramers-Kronig relations 
[6]. An expansion of the field operators is performed which is based on the 
Green function of the classical Maxwell equations and preserves the equal- 
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time canonical commutation relations. 

Suttorp and Wubs in the framework of the damped polarization model, 
have quantized the electromagnetic field in an absorptive medium with spa- 
tial dependence of its parameters [7]. They have solved the equations of 
motion of the dielectric polarization and the electromagnetic field by means 
of the Laplace transformation for both positive and negative times. The op- 
erators that diagonalize the Hamiltonian are found as linear combinations 
of canonical variables with coefficients depending on the electric suscepti- 
bility and the dielectric Green function. Also the time dependence of the 
electromagnetic field and the dielectric polarization are determined. 

The macroscopic description of a quantum damped harmonic oscillator 
with frequency u> is represented in terms of the Langevin equation [8,9,10]: 

x+ / dt'/i(t - t')x(t') - ul / dt'vit - t')x(t') = F N (t). (1) 

JO JO 

The coupling with the heat bath in the microscopic theory corresponds to 
two types of forces in the macroscopic equation of motion of a damped har- 
monic oscillator. The forces of the first type are obtained from some memory 
functions n and v. The second type force is the noise force F^it). These 
two types of forces have a fluctuation-dissipation connection and both are 
required for a consistent description of a dissipative quantum system. 

Matloob has quantized the macroscopic electromagnetic field in a linear 
isotropic permeable dielectric medium by quantizing the Langevin equation 
and associating a damped quantum harmonic oscillator with each mode of 
the radiation field [10]. There are some other approaches to quantizing the 
electromagnetic field and interested reader is referred to [11-21]. 

In this paper by modeling a polarizable and magnetizable medium with 
two massless quantum fields, namely " E and M quantum fields" , we come to a 
completely systematic and consistent method for quantizing electromagnetic 
field in such a medium. We propose a general method in which the medium 
is included in quantization process. In fact, the coupling of electromagnetic 
field with the medium in a microscopic level, is modeled macroscopically by 
replacing the medium with "E and M quantum fields ". The underlying 
mechanism of interaction between electromagnetic field and the medium in 
microscopic theory may be much more complicated than that implied by 
this model which simply replace the medium with a collection of harmonic 
oscillators. The present model is a generalization of Caldeira-Leggett model 
where for dissipative quantum systems, they model the environment by a 
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collection of harmonic oscillators [22,23]. In the Caldeira-Leggett model, the 
environment's Hamiltonian is 

^ = B^ + ^X], ( 2 ) 

where m n ,x n ,p n and u> n are mass, position, momentum and frequency of 
rath oscillator respectively. It should be noted that the properties of the 
environment may in some cases be determined on the basis of a microscopic 
model which is not necessarily the Caldeira-Leggett model. As an example we 
mention an Ohmic resistor which as a linear electric element should be well 
described by a Hamiltonian of the form (J2J). On the other hand the underlying 
mechanism leading to dissipation in a resistor may be more complicated than 
that implied by the model of a collection of harmonic oscillators. 

The main feature of the present approach is that "E and M quantum 
fields" are to describe the electric and magnetic properties of the medium 
macroscopically. It is clear (at least in the present work, based on the results), 
that this assumption is an effective one which can be compared with the 
phenomenological approaches to the problem. 

In the Green function method and schemes applied in references jSj, [TU] . 
[TB^ . jEj, [21] the noise current and polarization densities are due to in- 
teraction of electromagnetic field with the medium. Since the kind of this 
interaction is not defined, the explicit forms of the noise current and polariza- 
tion densities are not known. It is clear that the dissipative character of the 
medium depends on electric and magnetic susceptibilities. Since the presence 
of the noise quantum fields are necessary for quantum dissipative systems, 
there should be a dependence between the strengths of the noise densities 
and the strengths of electric and magnetic susceptibilities. Specially, when 
the medium becomes a non dissipative one, i.e., when the imaginary part of 
the Fourier transform of dielectric function tends to zero, the strengths of the 
noise densities should clearly tend to zero. In the Green function method and 
schemes applied in jH] JTD| , [IS] , [El , [21] , the relation between the strengths of 
the noise densities and the strengths of electric and magnetic susceptibilities 
is not clear. Therefore it is not clear how these noises tend to zero when the 
medium becomes a non dissipative medium. 

The constitutive equations of the medium which relate electric and mag- 
netic polarization densities to electric and magnetic fields respectively, should 
be treated as a consequence of interaction between electromagnetic field and 
the medium. But in the mentioned methods, the constitutive equations can 
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not be derived from Heisenberg equations of motion. In the present ap- 
proach, the coupling of electromagnetic field with the medium is known and 
explicit forms of noise densities are given. Therefore the relation between 
the strengths of noise densities and the strengths of electric and magnetic 
susceptibilities is clear. Also, in this approach, the constitutive equations 
can be derived from the Heisenberg equations of motion. The key role is 
played by what we have introduced as coupling functions, which couple the 
electromagnetic field to "E and M quantum fields". The electric and mag- 
netic susceptibilities are defined in terms of the coupling functions. Also, 
the noise polarization densities are described in terms of the coupling func- 
tions and the creation and annihilation operators of the "E and M quantum 
fields". The coupling functions are common factors in the noise densities 
and the electric and magnetic susceptibilities. It can be shown that for a non 
dissipative medium, the noise densities vanish as expected. 

In the damped polarization model the polarizability property of the di- 
electric, is described by a quantum field Y and the dissipative property of 
the medium is described by introducing a heat bath which is independent 
of the field Y. The heat bath is a continuum of Klein - Gordon fields with 
a continuous frequency range. The heat bath interacts with the medium in 
a suitable way. In this model, the magnetic property of the medium is not 
included [5,7]. 

In the present approach, the magnetic property of the medium is included 
and the polarizability and the dissipative properties of the medium are de- 
scribed only in terms of a single quantum field (E quantum field). Also, the 
electric susceptibility is defined in terms of a coupling function which couples 
the "E quantum field" to electromagnetic field. Finally, the electric polar- 
ization field of the medium is defined in terms of the coupling function and 
the creation and annihilation operators of the "E quantum field". 



2 Quantum dynamics 

Quantum electrodynamics in a linear polarizable and magnetizable medium 
can be accomplished by modeling the medium with two independent quantum 
fields which interact with the electromagnetic field. One of these quantum 
fields namely " E quantum field" , describes the polarizability character of the 
medium and interacts with the displacement field D through a minimal cou- 
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pling term. The other quantum field, namely " M quantum field" , describes 
magnetizability character of the medium and interacts with magnetic field 
through a dipole interaction term. The Heisenberg equations for eletromag- 
netic field (system) and " E and M quantum fields" (environment), lead to 
both Maxwell and constitutive equations. The constitutive equations relate 
the electric and magnetic polarization densities to the macroscopic electric 
and magnetic fields, respectively. 

The vector potential of the electromagnetic field in Coulomb gauge can 
be expanded in terms of the plane waves as 



where u>^ = c\k\ and Sq is the permitivity of the vacuum. The unit vectors 
e*£ A , (A = 1, 2) are polarization vectors and satisfy 

e^-k = 0. (4) 

These recent relations, guarantee that the vector potential satisfies the 
Coulomb gauge V • A(r, t) — 0. 

Operators ag A (£) and o>\At) are annihilation and creation operators of 
electromagnetic field and satisfy the following equal-time commutation rules 

[a^t),ai, x ,(t)] = Kk-k')5 xx , (5) 

The conjugate canonical momentum density of the electromagnetic field 
7rp(r, t) and also the displacement vector operator D(f,t) are by definition 
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(6) 

From this definition it is obvious that V • D — 0, which is the Gauss law 
in the absence of external charges. The commutation relations (JHJ lead to 
the commutation relations between the components of the vector potential 
A and the displacement vector operator D as 

[Atftfy-Dji?,*)] = [MrA^j(P^)]=^ij(f-P), (7) 
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where 5^{f—r') = ^3 / d 3 ke lk '^ r '- ) (^j — * s ^ e transverse delta func- 
tion with the following properties: 

1. Let us define the transverse and longitudinal components of an arbitrary 
vector field F(r,t), as 

F ± (r, t) = F(f, t) - f d 3 r'V ■ F(P, t)VG(f, P), (8) 
pW(r, t) = J d 3 r'V ■ F(P, t)VG(f, ?), (9) 

respectively. Wherein G(f, r') is the Green function 

G(r,P) = - 1 (10) 
47r|r — r I 

Then as the first property, one can easily show that 

F^f,t)=j2f d^f-P^P). (11) 
i=i 

2. The second property is the transversality of 5y (r — r') 

E^(r-r') = 0, J = 1,2,3. (12) 

The Hamiltonian of the electromagnetic field inside a magneto-dielectric 
medium, can be written as 

= E / ^ 3 ^4( £ ) a feA(*)> (13) 

where /xq is the magnetic permitivity of the vacuum and we have applied the 
normal ordering to operators &r A W an( l a k\(^)- 

Now as mentioned, we model the medium by two quantum fields, E 
and M, which describe the electric and magnetic properties of the medium 
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macroscopically. Therefore, the Hamiltonian of the medium can be written 

as 

H d = H e (t) + H m (t), 

3 

He{t) = ]T f d 3 q f d 3 khu$dl{k,q,t)d v {k,q,t), 

v=l J J 

H m {t) = Yl ! <Pk?w £ bUk,q,t)b v (k,#,t). 

u=l J J 

(14) 

where cug is the dispersion relation of the medium and H e and H m are the 
Hamiltonians of the " E and M quantum fields" , respectively. The operators 
d v {k,q,i), d\,{k,q,i), b u (k,q,t) and bl(k,q,t), are annihilation and creation 
operators of the " E and M quantum fields", respectively. We impose the 
following equal-time commutation relations on these operators 

[d v {k,q,t),dl{k',q',t)\ = 5 vvl 8{k-k')5{q-^), 
[b v $,q,t),bl<fi,4,t)\ = 8 vv ,8{k - k')S(q- q 1 ). (15) 

Now let us define the electric polarization density operator of the medium as 

P(r,t) = £ / -j=L= I d 3 k[f(^ k ,f)d u (k,q,t)e^+H.C]vM, 
u =i J J(2tt) 3J 

(16) 

where for convenience H.C stands for Hermit ian Conjugation and 
v v (q) = e v $, for v = 1,2, 



v 3 (q) = q = ||j , for v = 3. 



(17) 



In a similar manner we define the magnetic polarization density operator of 
the medium as 

M(r,t) =i£ I -£±= I d 3 k[g(^ k ,^b u (k,q,t)e^ - H.C.]s u (q), 



(27T) 



s v {q) = q x e V(i , for v = 1,2, 

Mo) = & for u = 3. (18) 
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In the polarization densities (fH)|) and ([18)1. the functions f(u>^, r) and g(iOj:, r) 
are the coupling functions between the electromagnetic field and " E and M 
quantum fields". The coupling functions are position dependent (indepen- 
dent) for a inhomogeneous (homogeneous) magneto-dielectrics. 

The polarization densities P and M, are defined on the basis of the fol- 
lowing physical assumptions: 

1- These densities should be Hermitian operators. 

2- For a linear magneto-dielectric medium, P and M, should be a linear 
combination of creation and annihilation operators of the medium. For a 
non- linear medium, P and M, may not have a linear expansion in terms of 
the creation and annihilation operators of the medium. 

3- Polarization densities should depend on the macroscopic properties of the 
medium. 

Macroscopic properties of the medium will be reflected in electric and 
magnetic susceptibilities. It is clear from relations (fTBj) and (fTHj) that the po- 
larization densities P and M depend on coupling functions f(u k ,r), g(u>k, r) 
and dispersion relation {uA of the medium. In the following we show that 
electric and magnetic susceptibilities are dependent on the coupling functions 
and dispersion relation. 

Now let us propose the total Hamiltonian as 

[D { r, I) - P(r, ,)]> , (V x AW, t) v x ^ . ^ 




2e Q 2/i 

(19) 

Using (J7J) and (jllj) . we can obtain the Heisenberg equations for A and D 

(20) 

^) = i hl H,5 ( r, t)] = VxV ;/>"'> - V x *V, t ), PD 

where P -1 and M -1 are transverse components of P and M respectively. The 
transverse electric field E- 1 , magnetic induction B and magnetic field H are 

a 7 d 
E ± = - — , B = V x 1, H= M, (22) 

at /i 
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so (|2(J|) and (J21)) can be rewritten as 

D = e E ± + P ± , (23) 

^ = VxH ± . (24) 

In absence of an external charge density, we have D" = £o-^" + -P" — 0, and 
we can define the longitudinal component of the electric field as i?" = — ^jL 
Combining (j2"U)l and ()21|1. we find 

1 <9 2 v4 

" V ^ + = ^°^T + ^oV x M\ (25) 

By using the commutation relations ()15|l. the Heisenberg equations for oper- 
ators <ij,(£;, g, t) and 6 y (fc, q, t) can be obtained as 

d v (k,q,t) = UH,d u {k,q,t)} 

= -iu % d v (k, q, t) + / / d 3 r'f*(^, P)E(P, t)e-^ ? ■ v u (q), 
HJ (27r) 3 j 

(26) 

and 

K(k,q,t) = -[H,b v (k,q,t)] 

= - iu % b v {k, q, t) + / / dV<7* (o^ P) e -^ P S(?, i) • s u (o), 

hJ (27r) 3 ^ 

(27) 

respectively. These equations have the following formal solutions 

d u (k,q,t) = d v (k,q,0)e* u i t + 

dt'e-^H-t') J d 3 r'f*{uj^?)e-^ p E{P,t')-v v {q), 

(28) 



^(2tt) 3 



' r * dt'e-^-V I d 3 r'g*(u) U , P)e-^ ? B(P, t') ■ s„(q). 



3 Jo 



(29) 
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By substituting (|28|l in (JTEJ), we find the polarization operator 

P(r, t) = P N (r, t) + e f *' dt'x e (r, |t| - t')^, ±0. ( 30 ) 

Jo 

where i£ = — ^ — ^- is the total electric field. The upper (lower) sign, 
corresponds to t > (t < 0) respectively. 
The memory function 

gyj- poo 

Xe{r,t) = — d\k\\k\ 2 \f{u h ^\ 2 smu % t t>0 
fe Jo 

Xe(r,t) = t<0 (31) 

is called the electric susceptibility of the magneto-dielectric which is defined 
in terms of dispersion relation uo^ and the coupling function f(us,r). The 
operator P N (f,t) in (|30j) is the noise electric polarization density 



pN(r,t) = J2 I "7== I rf 3 M/K,r1<(fc,g,0)e-^ t+ ^>//.a]^(g). 
v =i J \/(2tt) 3J 

(32) 



By using the equation (}3*T]) . we can obtain the following important relations 
in frequency domain 

'"fel-ll^fl/f")!- (33) 
te\zj?,u)) = r- d\k\\k\*\f{^?)\*^- 2 , (34) 

fc 

where 

X e (r"»= / cft Xe (r,i)e^ (35) 
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A feature of the present approach is its flexibility to choosing an appropriate 
dispersion relation such that Jj|r > 0. When a dispersion relation 

is given, then knowing the susceptibility, we can obtain the corresponding 
coupling function from ()33|) easily. The sign of the left hand side of the re- 
lation (}3*3*j) should be positive because it is the imaginary part of the electric 
susceptibility in the frequency domain and this imaginary part is also con- 
nected to losses in the medium, so it is necessarily positive, otherwise x (f, ui) 
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should be discarded, as it would be unphysical. Therefore for consistency one 
must choose a dispersion relation which is a strictly increasing function of 
\k\. Here we have chosen the simplest one, i.e. a linear dispersion relation. 
It is remarkable to note that for a given susceptibility x (r, lo), by using (JT5j) 

and the definition of Pjy in (JH2*j) . one can show that for any choice of 
and f(uj,r) satisfying (j3HJ), the following commutation relations 

fie 

[£^(r,u;),P^(rV)] = — ImW(f, u;)]<%<5(f - ?)8{u - u') (36) 

71" 

between the Fourier transforms of the components of the noise polarization 
density do not change. These commutation relations are identical with those 
in reference jT3| and lead to the correct commutation relations of electromag- 
netic field operators. So various choices of and f(u,r) restricted by 
(I33j) . do not affect the commutation relations (J36j) and also the commutation 
relations between electromagnetic field operators. This means that there are 
many models with the same electric and magnetic properties which can be 
taken as an environment with different dispersion relations. In this point 
of view w(|fc|) and f(u,r) are two free parameters of our model up to the 
relation (}33|) . As was mentioned, we choose the simplest dispersion relation 
k>(|fc|) = c\k\, where c, velocity of light, as the proportionality coefficient, is 
just for simplifying the calculations. It is clear from equation (}3*3*j) that other 
choices of the dispersion relation, just lead to a redefinition of the coupling 
function f(u,f) and also more difficult mathematical expressions. 

For the choice = c\k\, the definition of the electric susceptibility in 

()31|) . becomes 
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Xe(r,t) = 0, t<0, (37) 

For a definite Xe(^ t), which is zero for t < 0, we can obtain the corresponding 
coupling function /(wr, r), in terms of Xe(r, t), as 

\f{Uk,r)\ 2 = -^j2j Q dtxe{r,t)smuj % t, > 0, 

k 

\f(^,r)\ 2 = 0, Wjg = 0. (38) 
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Similarly, by substituting (|29|) in (|18|) . we obtain the following expression for 
the magnetic polarization density M(r, t), 

M(r, t) = M N {f, t) + — /" dt' Xm {r, \t\ - t')B(r, ±t'), (39) 
where Xm is the magnetic susceptibility of the magneto-dielectric 

,_ . 87ryU, f°° , 21/ 

Xm(r,t) = J o du^\g{u^r)\ sin^t, t > 0, 

X m (f,t) = 0, £<0. (40) 

If we are given a definite x m (r, £), which is zero for t < 0, then we can obtain 
the corresponding coupling function g{uj^r) in terms of X m (r, £) as 

l^fc^)! 2 = 4 7r 2 //QU; 2 y dt Xm(r,t) sin wgt, > 0, 

|p(wff,»0r = 0, w ? =0. (41) 

The operator Mjv(r, £), is the noise magnetic polarization density 

Miv(r,t) = 
i/=l J J w(27t) 3 

(42) 

It is remarkable to note that the constitutive equations ("2HJ), (|3~0j) and 
together with Maxwell equations are obtained directly from the Heisenberg 
equations applied to the electromagnetic field and the quantum fields E and 
M. The explicit forms of the noise polarization densities are given by (p"""") 
and (|42)1. The coupling functions / and g, are common factors in the noise 
densities Pn and Mn and the susceptibilities Xe, Xm- So it is clear that the 
strengths of the noise fields are dependent on the strengths of the suscep- 
tibilities Xe and x m , which describe the dissipative character of a magneto- 
dielectric medium. 

For a homogeneous medium, the coupling functions f{oJ%, r) and g(ui£,r), 
are position independent. In this case, from ([3*7)1 and (jUJJ), we deduce that 
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Xe and Xmi are a l so position independent. Substituting (JHUj) and (|39p in the 
wave equation ()25|). we find 

V x y dt' Xm (\t\-t')V x A(r,±0 = -Qf(r,t) + ^ V x M^(f,t), 

(43) 

where c 2 = and the upper (lower) sign corresponds to t > (t < 0) 
respectively. 

The equation (}4*3*|) is the Langevin-Schrodinger equation [8] for vector the 
potential A, wherein, the explicit form of the noise current density is known. 

The quantum Langevin equation can be considered as the basis of the 
macroscopic description of a quantum particle coupled to an environment or 
a heat bath |B],|S]. 



3 Solution of the Heisenberg equations 



In this section, we solve the Heisenberg equations for the vector potential. 
Let us denote the Fourier transform of the vector potential A(r, t) by A(q, t), 
so 



A(r,t) 



(2*0* 



d 3 qA{q,t)e^, 



(44) 



from (PJ, it is clear that 



^ ^ 2e w f 

The wave equation (|43[) can be written in terms of A(g, t) as 



(45) 



A + ^ A± aiL 

d 3 k 



— £ 

£ A =l 



'2vr) 



^(1*1 - oi(9, ±0 - w ? y ^'xmdti - f ±0 

[^-/(^(fc, 0)e-^e^ - WjE /*(ufc)4(£, -9, 0)e^ 
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where u$= c\q\. 

This equation can be solved using the Laplace transformation method. 
For any time dependent operator g(t), the forward and backward Laplace 
transformations are by definition 

poo 

g f ( 8 ) = / dtg(t)e- st , (47) 
J o 

alld ,00 

g b (s) = / dtg(-t)e- s \ (48) 



respectively. Let Xe(s) and Xm(s) be the Laplace transformations of Xeif) 
and Xm(t) respectively. Then Af(q,s) and ^(q, s), i.e., the forward and 
backward Laplace transformation of A(q,t), can be obtained in terms of 
Xe(s) and Xm{s) as 

1 U 
S 2 S 2 Xe(s) -u}Xm(s)~ ^ 

£ J \ S ± UJ S ) [s 2 + Ul + S 2 US) - u}Xm{s)\ ^ 

(s =p kjj) [s 2 + u;| + s 2 Xe(s) - ujjxm(s)] qX 
V- f J3u( 9(^k)h(k,q,0) 



\U 9 fey fci ( S ±^)[ S 2 + w |+ s ^ e ( S )- w |x m ( 5 )] 9A 

+ ( STi ^)[ s 2 + ^ + s ^ e ( s )_ a; ^ m(s) ] -^-r i J 

where the upper (lower) sign corresponds to Af(q,s) (^(g, s)) respectively. 
Now taking the inverse Laplace transformation of Af(q,s) and A 6 (g, s), we 



(46) 
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obtain a complete solution for A(f, t) 

2 



A(r,t) = E / d% 9 , 9 ^ 3g , [Z±(^t)e^% x (0) + H.C.]e q ^ 

A=l \ £ W 5 

£ o x =i J tJ(2tt) 3J 

(50) 



where the upper (lower) sign corresponds to t > 0(t < 0). The functions 

Z + (UJ ? , t), Zjpft -t), C+( W fc> U)ft t), C-(^fc, ^q, ~t), Uq*, t), 7/_(Wjf, LO S , -t) 

are given by 



Z±(u& ±t) = L 



-i 



[s + sxe{s) =F 



>7±(<*>£,W ?J ±t) = #(^fc)£ 



-1 



(S ± KJj) [S 2 + w| + S 2 Xe(s) - w|Xm(s)] 
(S ± ^)[S 2 + + S 2 Xe(s) ~ ujXm(s)] 



(51) 



for t > and L _1 /(s) is the inverse Laplace transformation of function /(s). 
The transverse component of the electric field can be obtained from E ± = 
— Having the vector potential A, the transverse component of electric 
polarization density (P ) and also magnetic polarization density (M), can 
be obtained easily from relations (JSOj) and (jHHJ). 

Taking the Laplace transformation of the constitutive equation the 
longitudinal component of the electric field can be written as 

E\r,t) = - — 



d 3 q 



to J J(2n) 



d*k[Q ± {u h t)f(u> n )<k(k, q, 0)e^+ #.C.]g, 

(52) 
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where Q + (u^,t), Q_(co>£, — t) are given by 



Q ^ ±t) = L ~ l {jrrm: 



»( S ^r (53) 

for t > 0. In the following we consider some important examples. 
Example 1: 

Let f{oJ%) = g(&j;) = 0, then from (p?T|) and (}4T)j) . we have Xe(t) = Xm(0 = 0, 
and from (joTj) . we find 

*) = Z.(u € , t) = e~^\ V± = C± = 0, (54) 

therefore in this limiting case, quantization of electromagnetic field reduce 
to the usual quantization in the vacuum as expected. 

Example 2 : 

Take Xeif) and Xm(t) as follows 



Xe(t) 



% < t < A, 
otherwise, 

< t < A 



Xm (t) = { (xS^+i)A 

[ otherwise, 

(55) 

where x° e i an d A are some positive constants, using (f3SJ) and f|41|) . we 
find the corresponding coupling functions as 

ftc 3 £ Xe Sin2 ~2~~ 



2 



,„ / x.2 _ ^ C Xm 8111 "IT 



(56) 



4ttVo(^ + 1)^ ^ 



and from (|50|) and we have 

P(r, t) = P^f, t) + t A dt'E(r, ±0, 

A J\t\-A 

M(r,t) = M N (f,t) + *° m dt'B(f,±t'), (57) 
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where Pjv(r, t) and M^{r, t) are the noise polarization densities and 
corresponding to the coupling functions obtained in (|56|) . 

In the limit A — > 0, the coupling functions ()56|) and the noise polarization 
densities, tend to zero, and the relations (|57jl are reduced to 



P(r,t) = e X° e E(r,t), 



M{^t) = X ™ B(r,t). 



(58) 



In this limit, the electric field and the polarization densities are purely trans- 
verse and 



Z-(u)fa t) = Z + (UJ$, t) = COSCUq-t — i, 



1 + y° 

Sill UJgt, 



r)±(^t) = C±K-,i) = 0. (59) 



Also the electromagnetic energy inside the dielectric is 

' 7 A(£0)-i'(£0)] (60) 



2(1 + 

where D is the Fourier transform of the displacement field. The energy 
given by l6T?j) is a constant of motion contrary to the vacuum expression 
JI^EqE 2 + ^}d 3 r, which is not clearly a constant of motion. This example 
shows that this model can be applied to a nondispersive magneto-dielectric 
medium. 

Example 3: Let Xeif) = (3u{t) and Xm(t) — 0, where u(t) is the step 
function 

' 1 t>0 
t < 

(61) 



u{t) 
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and j3 is a positive constant, then using (|38|) and (|41j) . we find 

2 hc 3 e /3 



Ak 2 ujI 

k 



(62) 



and accordingly we can rewrite (j4T)j) as 



— I 



hc 3 p 



d 3 k 



(63) 



which has a dissipative term proportional to the first time derivative of dy- 
namical variable A. From (JoTj) . one can obtain 

Z±(u>$, t) = e T ^*[±-^- sin Vt^t + cos fi^t — — p sin fi^t], 



\ 



hc 3 fte 
An 2 , 



T— ^ 



T g f (-f + tfy)e ±tf V (f + zggje^V , 



(64) 



where fig- 
is 



w| — ■ The asymptotic solution of A(r, t) in large-time limit 



d 3 q 



fic 3 (5 * 



d x (k,q,0)e- 1 ^ 

<4 



-B.C. 



(65) 



Using 1)5 2j) . it is easy to show that the longitudinal component of the electric 
field in the limit t — > ±oo is 



#11 (ft) = 



hc 3 /3 r d 3 q r d 3 k 



An 2 En 



(2tt) 3 



/LUt 



d 3 (k,q,0) 
P =F i^n 



(66) 
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Example 4: A simple model for x e (s) 

If we neglect the difference between local and macroscopic electric field for 
substances with a low density, then the classical equation of a bound atomic 
electron in an external electric field is 

f + 7^ + uj 2 f= -—E(t), (67) 

where the influence of the magnetic force has been neglected compared to 
the electric force. The parameter 7 is a damping coefficient and the force 
exerted on the electron due to atom is taken to be simply a spring force with 

frequency loq. If E(s) and r(s) are the Laplace transformations of E(t) and 
r{t) respectively, then from (jBTf) we can find 

r(s) 



S z + 7s + UJq 

Now let there be iV molecules per unit volume with z electrons per molecule 
such that fj electrons of any molecule have a bound frequency Uj and a 
damping coefficient jj. The Laplace transformation of the polarization den- 
sity is 

~-> Ne 2 f- ~u 

P(s) = — E ~2 2 E ( S )- ( 69 ) 

If Uj and 7j are identical for all of electrons, then from (JfiUj) we can write 

v (,) - ^ u 2 - — 
Xe{s) — ; : : ; — -, Co> p - 



s + 7s + ujq v me 



, 2 



9 7* Sin fnt , . 9 9 7 

Xe (t)=^e-- °-u(t), u 2 =u 2 --L, (70) 

where is the step function defined in (J6l"|) . We can obtain the coupling 
function f(u%) from ([HEj) as 



2 



7 7 



^ [ ^ + (i/ - ^) 2 1 + (^0 + ^) 2 

If 7 = 0, then the dielectric substance is a nondissipative one and the coupling 
function takes the form 

l/K)| 2 = ^^^o-^). (72) 
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In this case, the noise electric polarization density is nonzero only for the 
resonant frequency (u = Uq) of the equation r + UqV = —^E e~ luJot . In 
the resonant case, energy of electromagnetic field will be absorbed by the 
medium. By using (JoTj) . we can obtain Z + (u$,t) and Z^(cu^, —t), as 

z±K-i ±() = L " { * + IZi +V+% + + %i ) ' (73) 

for t > 0. The equation (J73j) can be solved by calculating the residues of the 
function 

[(sTt^)(s 2 + u 2 ) + su 2 p }e ±st 

by extending s to the domain of complex variables. Similarly C± can be 
obtained from (J51|) with f{uj£) given by (|72j) . In this case we have 



Q±( U * t ) = . .2 , .2 fc ..2 e 



iujrt 

K 




The longitudinal component of the electric field can be obtained from (|52|) . 
using and Q±(^,t) given by (J72j) and (J7HJ) respectively. 

If 7 7^ 0, the substance is of dissipative kind and A{r, t), can be obtained 
from ()5(jj) with /(o;^) given by (|7T|) and Xe( s ) given by (|7UJ) . In this case the 
longitudinal component of the electric field can be obtained from (|52|). This 
example shows that this model of quantization of the electromagnetic field 
is applicable to both dissipative and nondissipative dielectrics. 



4 Concluding remarks 

By modeling a linear, polarizable and magnetizable medium with two quan- 
tum fields E and M, electromagnetic field is quantized in the medium, consis- 
tently and systematically. There are many models for the environment with 
the same electric and magnetic properties. In other words, there are many 
dispersion relations and coupling functions leading to the same suceptibili- 
ties of the medium. Therefore one can take the simplest physical dispersion 
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relation. Once a dispersion relation is defined, then for any definite magneto- 
dielectric medium, i.e., Xe{t) an d Xm(t) are known functions, one can find the 
corresponding coupling functions fioj^r) and g(uj:,r). The coupling func- 
tions describe the electric and magnetic properties of the medium macro- 
scopically. These functions also couple electromagnetic field to the medium 
through the quantum fields E and M. The explicit forms of the noise den- 
sities are derived. Since the coupling functions are common factors in the 
noise densities and susceptibilities, the relation between the strengths of the 
noise densities and susceptibilities is clear. In this approach, both Maxwell 
and constitutive equations are obtained as Heisenberg equations of motion. 
In the limiting case, i.e., when there is no medium, the approach tends to 
the usual quantization of the electromagnetic field in vacuum as expected. 
This model of quantization is applicable to both dispersive and nondispersive 
magneto-dielectrics. 
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